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ABSTRACT 

Accurate rotational-vibrational line lists for and in their ground 

electronic states are computed. The line lists are produced using an empirical potential 
energy curve obtained by fitting to the experimental transition frequencies available 
in the literature in conjunction with an accurate, high level ab initio dipole moment 
curve. In these calculations the programs DPotFit and LEVEL 8.0 were used. The 
new line lists reproduce the experimental wavenumbers with a root-mean-square error 
of 0.004 cm“^. The line lists cover the frequency range 0-51000 cm“^, contain almost 
700 000 lines each and extend up to a maximum vibrational level of v=66 and a 
maximum rotational level of J=357. They should be applicable for a large range of 
temperature up to, at least, 5000 K. These new line lists are used to simulate spectra 
for PN at a range of temperatures and are deposited in the Strasbourg data centre. 
This work is performed as part of the ExoMol project. 

molecular data; opacity; astronomical data bases: miscellaneous; planets and satellites: atmospheres; stars: low-mass 


1 INTRODUCTION 


Phosphorus nitride, PN, is a stable, strongly bound, cl osed shell diatomic molecule. Its electronic ground state has 
symmetry and its dissociation energy is 51940±170 cm“^ ( Gingeri^l 19691 : Cazzoli et al.l[20oJ ). It has no low-lying electronically 
excited states: the lowest excited state is the a^E"*" state at about 25000 cm~^ (for near-equilibrium geometries, ~1.5 A) 


while the lowest-lying singlet is the A ^11 state at about 40000 cm ^ ( Grein fc Kaour IPS,*! : de Brouckere et Tl 199 J) . 


PN is an astrophysi cally importa nt molecule used to probe different regions of the interstellar medium (ISM). PN was first 
observed in the ISM by Ziurvsl (Il987 ). where it was noted that PN was more abundant than o riginally predicted. It was also 


observed in the warm star forming regions Ori (KL), W51M, Sgr B2 by Turner fc Ballvl (119871 ) as well as in warm molecular 
clouds of the star forming regions M17SW and DR 120H by Turner et al.l ( 199fll ). sug gesting high temper ature chemistry is 
important for its formation. PN was detected in the c arbon-rich source CRL 2688 by Milam et al. ( 20081 ). Most recently it 
was observed i n the Lynds 1157 B1 sho cked region bv lYamaguchi et al.l (|201ll ) and in the outflow from the pro tostar IRAS 


20386 -h6751 bv lYamaguchi et al.l (l2012l) . The chemical mechanism of the PN formation was recently studied bv IViana et al 


1 2 OO 9 I) . As yet PN has not been observed in hotter bodies; however Visscher et al. 1 20061) studied phosphorous and found that 


in hotter objects, such as brown dwarfs, P-bearing gases such as PN become increasingly important at higher temperatures. 
Similar results can be expected from exoplanets. 


were unaeriaKen in ine suD-miiiimetre, raaiowave, microwave ana mrra-rea irequ 
Hoeft et al.l 1972: Wvse et al. 1972: Coauart & Prudhommel 1980l. 1981: Maki & 

encv ranges im.avmonaa v ivie nDerei 
Lovadll981: Ghosh et al. 198ll: Verma et al. 

19871: Le Floch et al. 1996: Cazzoli et al. 20061). Some molecular and spectroscopic constants for PN computed ab initio were 

reported bv Ravmonda & Klemperer I1971^: Wvse et al. 1 

1972h: Ghosh et al. ( 

198lh: Saraswathv & Krishnamurtv (1987): 

Ahmad & Hamilton (Il995). Several attempts have been made to construct potential energy curves of PN JMaki & LovasI 

198ll: Grein & Kaouilll983: Wong & Radorn 199fil: McLean et al. 19921: de Brouckere et al. 1993: Kemenv et al. 

2003). The 


molecular beam electric resonance measurements. 

The aim of this work is to produce accurate rotational-vibrational (ro-vibrational) catalogues of transitions (so called line 
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2 Yorke et al 


lists) for the two stable isotopologues of PN, and in their ground electronic states. To this end an ab initio 

dipole moment curve (DMC) is used in conjunction with an empirically rehned potential energy curve (PEC). Our line lists 
should be complete enough to describe the absorption of PN up to 5000 K. This work is performed as part the ExoMol project, 
which aims to provide line lists for all the molecular transitions of importance in the atmospheres of (exo-)planets and cool 


stars I Tennyson fc Yurchenko 


BeH MgH, CaH dVadin et al 
2014) . 


201^. The ExoMol methodology has been already applied to a number of diatomic molecules: 


lop 

20121 1 ■ SiO (iBarton et al.ll2013 l. NaCl and KCl (iBarton et al.ll2014l l. and AlO llPatrascu et al. 


2 SOLVING THE NUCLEAR MOTION PROBLEM 

A spectroscopic line list comprises transition frequencies (i.e., line posit ions), Einstein coefficients as well as energies, degen¬ 
eracies and quantum numbers of the states involved in the transition ( Tennyson fc Yurchenkol 2012 1. When combined with 
a partition function, a line list can be used to produce spectra at a given temperature, both in absorption and emission. 
Each record in the line list corresponds to a transition between two energy levels, in this case ro-vibrational. These levels 
correspond to characteristic nuclear motions of the molecule, and thus to compute a line list one has to solve the corresponding 
Schrodinger equation and obtain the energy levels required. 

The general form of the radial Schrodinger equation for a diatomic molecule with atoms A and R in a electronic 
state is given by 


2fj, dr 2 


V{r) + 


2fj,r' 


:J{J + 1) 


‘4>v,j{r) = Ev,j'ipv,j{r), 


( 1 ) 


where V{r) is the potential energy curve, and il!v,j are, respectively, the bound state eigenvalues and eigenvectors and 
fj. is the reduced mass given by 

mAniB 




niA + ms 


In this work for all nuclear motion calculation we used the following (atomic) masses (IPfeiffer et al 
30.97376200041 u, m(^'‘N) = 14.00307400461 u and m(^®N) = 15.00010889818 u. 


2012h : 


Einstein A coefficients in s 


for rovibrational transitions in E electronic terms are given by I BernathI 2005h 
647rH0-3® S{J',J") _3, 


A = 


■ V \iYv' 


( 2 ) 

m(®ip) = 


(3) 


3h 2J' + 1 

where J' and J” are, respectively, the upper and lower rotational angular momentum quantum numbers, h = 6.62606957 x 
10“^^ erg s is Planck’s constant in cgs units, S{J', J") = max( J', J") is the Honl-London rotational intensity factor, v is the 
transition wavenumber in cm“^ and M(r) is the dipole moment curve (DMC) in debyes. 

Integrated absorbtion cross-sections I (also referred to as integrated absorption coefficients of simply as “line intensities”) 
in cm molecule”^ can be obtained by 

1 5rns(2J' + 1) g-£"/('=BT) ^-hcO/(kBT)^ 


I = 


1 — e 


(4) 


SttcO^ Q{T) 

where c = 2.99792458 x 10^° cm s“^ is the speed of light in cgs units, E” is the lower state energy, fcs is the Boltzmann 
constant, T is the temperature, Q{T) is the partition functions, and Qns is a nuclear spin factor. 

To solve the radial Schrodinger equation © and co mpute the ener gy levels, spectroscopic constants and Einstein coeffi¬ 
cients required we make use of the program LEVEL 8.0 IILe Rovll2007l '). The PEC and DMC are two prerequisites for these 
calculations and can be obtained from first principles by solving the Schrodinger equation for the motion of the electrons. Since 
ab initio PECs are not usually capable of providing accuracies better than a few cm~^, they are commonly refined by fitting to 
experimental energies or transition frequencies, subject to availabi lity. Conversely, the accuracy of ab initio DMC can be hig h 
enough to be competitive with , or even bett er than, experiment ( Lvnas-Grav et al. 19951 : Lodi et al. 2011 : Tennyson 2014). 
We use the DPotFit program ( Le Rov 2006ll to empirically refine the ab initio PEC by fitting to experimental frequencies 
collected from different literature sources, see Table [T] 


3 POTENTIAL ENERGY CURVE 

We produced in this work three potential curves for PN: an ab initio one named PEC-A; a semiempirical one named PEC-S 
which modihes the ah initio PEC with a single, physically-motivated, empirical parameter; and, finally, an empirical PEC-R 
obtained by htting the available experimental data. The three PECs are discussed in the following sections. 
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Exomol molecular line lists VI: PN 3 


Table 1. Summary of experimental spectroscopic data available for PN used in our fits. 


Ref. 

Isotopes considered 

max V 

max J 

Wavenumber range cm ^ 

Number of lines 

Cazzoli et al. (2006) 

31p 14[^ 

0 

17 

3.1-26.6 

24 

Cazzoli et al. (2006) 

31pl5N 

0 

15 

3.0-20.9 

9 

Maki &; Lovas fl981) 

31p 14]^ 

4, Av = 1 

53 

1217.22-1318.60 

22 

Wvse et al. ('1972) 

31p 14]\^ 

4, Av = 0 

8 

3.1-12.5 

19 

Ahmad Hamilton ('1995') 

31p 14j,j 

1, Aj; = 1 

33 

1273.30-1368.49 

62 

Ghosh et al. ('1981')“ 

31p 14j,j 

11 

0 

1323.17 - 13786.95*' 

11 


“ Vibrational band centres derived from the hi gh resolution study of ^n-^S+ transition. 
^ Range for the lower ( S+) state term values lIGhosh et alJll981^ . 


Table 2. Vibrational (J = 0) energy term values obtained with different electronic stru ctur e methods and basis sets. The second column, 
labelled ‘obs.’, reports experimental (observed) values from I Ahmad &: HamiltonI lll995l l and llrikural ||2007) 1. other columns show observed 
minus calculated. The stated uncertainty of experimental levels is 0.001 cm“^ for v = 0—4 and 0.05 cm“^ for v = 5-11. The energy curve 
corresponding to the results in the last column is referred to as PEC-A in the text. In all cases a relativistic correction curve computed 
with the MVDl Hamiltonian at the MRCI/aug-cc-pV6Z level was included. See Section Ifi.ll for detailed explanations. All values are In 
cm“^. 


V 

obs. 




obs. - calc. 






awc5z 

awc5z 

awc5z 

awc5z 

awc5z 

ac6z 

ac[56]z 



CCSD 

CCSD(T) 

CASSCF 

MRCI 

MRCI+Q 

MRCI+Q 

MRCI+Q 

0 

666.79 

-40.28 

-6.10 

9.72 

-5.85 

-1.77 

-2.55 

-3.37 

1 

1989.94 

-122.06 

-18.78 

28.93 

-17.67 

-5.37 

-7.70 

-10.16 

2 

3299.24 

-205.74 

-32.15 

47.80 

-29.71 

-9.07 

-12.97 

-17.07 

3 

4594.64 

-291.40 

-46.24 

66.34 

-41.97 

-12.88 

-18.35 

-24.09 

4 

5876.14 

-379.07 

-61.11 

84.57 

-54.44 

-16.80 

-23.84 

-31.24 

5 

7143.83 

-468.69 

-76.61 

102.66 

-66.98 

-20.68 

-29.31 

-38.38 

6 

8397.42 

-560.59 

-93.07 

120.35 

-79.86 

-24.79 

-35.02 

-45.77 

7 

9637.02 

-654.68 

-110.35 

137.83 

-92.93 

-28.97 

-40.82 

-53.27 

8 

10 862.51 

-751.11 

-128.59 

155.04 

-106.26 

-33.31 

-46.79 

-60.96 

9 

12 073.43 

-850.38 

-148.20 

171.60 

-120.27 

-38.22 

-53.36 

-69.26 

10 

13 270.68 

-951.60 

-168.27 

188.49 

-134.00 

-42.74 

-59.56 

-77.22 

11 

14 453.74 

-1055.34 

-189.33 

205.30 

-147.92 

-47.35 

-65.86 

-85.31 


3.1 The ab initio PEC-A 


We computed PECs using the well-known CCSD (coupled cl uster single and double e xcitations) and CCSD(T) (coupled cluster 
single, double and perturbative triple excitations) methods ( Bartlett fc MusialboOT ) as well as with CASSCF (complete active 
space self consistent field), MRCI (i nternally contracte d multi-reference configuration interaction) and MRCI-I-Q (MRCI with 
renormalized Davidson correctio n) ( Szalav et al. 2012h : the -l-Q correction was computed using the so-called ‘fixed’ reference 
function, see Szalav et al.l ( 20121 ) and the MOLPRO manual for an explanation. CASSCF and MRCI calculations used the full 
valence reference space and all electrons ap a rt from the phosphorus Is were correlated. We used cor relation-consistent basis 
sets of the Dunning family ( Dunning 19891 : Woon fc Dunning Jr. 19951 : Peterson fc Dunninel 2002 ) and in particular aug- 
cc-pwCV5Z and aug-cc-pC6Z, for which we use the shorthand notation aw c5z and ac6z. We extrap olated the energies using 
the formula = E^o + A{n + |)~^, which is known to perform very well ( Feller &: Petersonll2013l ). A relativistic correction 
curve was computed as expectation value of the mass-velocity one-electron Darwin operator (MVDl) with the MRCI/ac6z 
electronic wave function and was always included in the calculations. Energies were computed from r = 2 .10 ap to r = 6.10 a p 
(1.11 to 3.23 A) in steps of 0.05 ao and from r = 6.5 ao to r = 16 ao is steps of 0.5 ao using MOLPRO ( Werner et al. 20121 ). 
Calculation of one MRCI energy in the largest ac6z basis set took about 25 GB of disk and 3.5 hours on a single CPU. 

Some indicative results are reported in Table (2] From these it is clear that errors in ab initio energy levels are very large. 
Let us first discuss the results in the awc5z basis set. CCSD performs by far the poorest of the methods considered; the (T) 
correction reduces CCSD errors by a factor about 6 but the errors are still very large. This seems to indicate that at the very 
least CCSDTQ (coupled cluster with full triple and quadruple excitations) should be used for sub-cm“^ accuracy, which is 
computationally unfeasible. CASSCF energy levels are slightly worse than CCSD(T) ones. MRCI energy levels are slightly 
better than CCSD(T) ones, with errors reduced by a factor 1.1 or so. Finally, MRCI-I-Q improves considerably upon MRCI, 
reducing errors by a factor about 3. They are still, however, considerably off, with absolute errors in cm~^ progressively 
increasing with the vibrational quantum number as about 4.15 x (u -|- 1/2). 
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4 Yorke et al 


Note that all methods, apart from CASSCF (which anyway is not expected to produce quantitative accuracy), produce 
energy levels which are too high, suggesting the outer wall of the potential curve is rising too steeply. Increasing the basis set 
to ac6z leads to higher energies levels and therefore to even larger errors; basis set extrapolation inevitably leads to errors 
which are larger still. 

We estimate the residual basis set incompleteness error in extrapolated ac[56]z energy levels as half of the difference 
between ac[56]z and unextrapolated ac6z value; this estimate is well approximated by the expression 0.84x (u + 1/2), indicating 
that residual basis set errors are about hve times smaller than the observed errors. The relativistic MVDl correction curve 
(included in all calculations reported in Table [2j improves agreement with experiment; its effect on vibrational energy levels 
V = 0-11 can be expressed to within 0.04 cm“^ by —1.38 x (n + 1/2). Further corrections not considered in this work such as 
higher relativistic effects due to Gaunt, Breit or quantum electrodynamics corrections, as well as adiabatic and non-adiabatic 
effects, are certainly considerably smaller than the included MVDl curve and therefore their inclusion at this stage would 
not lead to an improvement of the observed residuals. We reach therefore the conclusion that the observed large errors in our 
best ab initio results are mostly due to incomplete treatment of electron correlation at the CCSD(T) or MRCI+Q level. 

Our highest-level ab initio PEC is the MRCI-|-Q/ac[56]z/MVDl, based on extrapolation of the awc5z and ac6z basis sets. 
This latter potential energy curve will be referred to as PEC-A and leads to a potential well depth De = 51 605 cm“^, which 
is lower by 335 cm“^ than the (rather inaccurate) experimental value 51 940T170 cm“^. 

It is worth noting that, almost certainly because of partial cancellation of errors, best agreement with experiment is 
obtained ab initio with the MRCI-|-Q/acw5z/MVDl PEC, at least for the range characterised by the available v = 0-11 
experimental data, which are sensitive to the PEC in the range r « 1.28-1.80 A. 


3.2 The semi-empirical PEC-S 


We considered a simple one-parameter semi-empirical modification of PEC-A as a way of improving its quality in a physically 
motivated way and with a view to provide a reasonable behaviour in the high energy range, experimentally unexplored. In 
particular, we investigated simple scalings of the Davidson correction contribution to the energies. The ah initio MRCI-I-Q 
energies can be broken down into a sum of three contributions: Emrci+q = Ecasscf + Ecorr + E+q, where Ecasscf is the 
energy of the reference CASSCF wave function, Ecorr = Emrci — Ecasscf is the MRCI correlation energy and E+q is the 
renormalized Davidson correction, which is given by E+p = Ec;orr(l — Co)/co where co is the expansion coefficient of the 
CASSCF wave function in the MRCI one ( Szalav et al. 2 OI 2 II . The Davidson correction tries to approximate the effect of 
quadruple excitations neglected in the MRCI calculation; in our calculation the cq expansion coefficient as a function of r 
has the value 0.966 for the near-equilibrium geometry r = 1.50 A, then decreases linearly up to about r — 2.0 A, assumes 
a minimum value 0.963 at r = 2.30 A and then slowly grows up to a dissociation value of 0.965. The small variation of 
Co with geometry and the fact co ~ 1 show that the CASSCF wave function is indeed a good reference function for all 
geometries. By contrast the expansion coefficient of the Hartree-Fock reference in the MRCI expansions is 0.91 at r = 1.50 A 
and then decreses monotonically, is 0.5 at r = 2.56 A and has an asymptotic value of about 0.28. Eventually we chose a 
simple direct scaling of the Davidson correction part of PEC-A by a factor a = 1.835. The optimal value of a was found by 
fitting the experimental J = 0 vibrational levels for v = 0-11 reported in Table [5] Equivalently, PEC-S can be written as 
Vs(r) = VaY) + 0.835 E+Q(r), where VA{r) refers to PEC-A. 

This one-parameter semi-empirical PEC, referred to as PEC-S below, reduces very considerably the residuals with ex¬ 
periment of PEC-A and leads to agreements of better than 0.30 cm“^ for to u = 0-8, while last three levels v = 9-11 show 
larger residuals of the order of 1.5 cm“^, see Table |3l The PEC-S leads to a potential well depth De = 51 742 cm“^, which is 
lower by 198 cm“^ but compatible with the experimental value, 51 940T170 cm“^. PEC-S was used as the starting point for 
our subsequent empirical fit. 


3.3 The empirical PEC-R 


We considered a fully em pirical potential energy function based on the Extended Morse Oscillator (EMO) function (ILe Rov 
2007 : Surkus et al. 1984 1 as given by 

Nb 


V{r) = De\l — e where /3(f) = 


rP _j_ 


(5) 


The parameters to be fitted are the equilibrium distance re, the potential well depth De, the coefficients pi and, finally, the 
exponent p. Note that one should constrain p > 0, /3i > 0 in order to insure l/(-|-oo) — De- In our fits we truncated the 

expansion of P{r) to t hree terms (N b = 2) for r < Ve- 

We used DPotFit ( Le Rovll2006l l to refine the expansion parameters by fitting to the experimenta l frequencies summa rised 
in Table[T] The main par t of the experimental data s et includes pure rotational transitions recorded by Cazzoli et al. ( 2006l l and 
infrared transitions from Ahmad fc Hamilton 1 1995l l which cover u ^ 4. In order to improve the sampling for higher vibrational 
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Table 3. Differences of experimental (‘obs.’) ro-vibrational energy levels with calculated ones (‘calc.’) for three potential energy curves 
X{r), X =A, S or R. PEC-A is an ab initio curve based on MRCl+Q/ac[56]z/MVDl; PEC-S is a one-parameter modification of PEC-A; 
PEC-R is an 8 -parameter empirical curve defined by Eq. l[5]l and table|4] Experi mental energy l evels are computed using the spectroscopic 
parameters by I Ahmad &: HamiltonI lll995h and include the zero point energy bv llrikur^ ll2007f ). 


V 

obs. 

obs. 

— calc. 


obs. 

obs. 

- calc. 


obs. 

obs. 

— calc. 




A 

s 

R 


A 

s 

R 


A 

s 

R 



o 

II 




J = 10 




J = 20 



0 

666.79 

-3.37 

0.02 

0.00 

752.99 

-3.53 

0.03 

0.00 

995.76 

-3.99 

0.04 

0.00 

1 

1989.94 

-10.16 

0.07 

0.00 

2075.53 

-10.33 

0.08 

0.00 

2316.58 

-10.80 

0.09 

0.00 

2 

3299.24 

-17.07 

0.10 

0.01 

3384.21 

-17.24 

0.11 

0.00 

3623.54 

-17.72 

0.12 

0.00 

3 

4594.64 

-24.09 

0.10 

0.00 

4679.00 

-24.27 

0.10 

0.00 

4916.60 

-24.76 

0.12 

0.00 

4 

5876.14 

-31.24 

0.06 

-0.01 

5959.88 

-31.42 

0.06 

-0.01 

6195.74 

-31.93 

0.08 

0.00 

5 

7143.83 

-38.38 

0.13 

0.14 

7226.93 

-38.58 

0.12 

0.12 

7461.00 

-39.16 

0.07 

0.07 

6 

8397.42 

-45.77 

0.04 

0.18 

8479.94 

-45.94 

0.06 

0.20 

8712.36 

-46.43 

0.12 

0.27 

7 

9637.02 

-53.27 

-0.07 

0.28 

9718.90 

-53.47 ■ 

-0.08 

0.28 

9949.49 

-54.04 ■ 

-0.09 

0.28 

8 

10862.51 

-60.96 

-0.29 

0.38 

10943.77 

-61.16 ■ 

-0.28 

0.39 

11172.63 

-61.71 

-0.27 

0.43 

9 

12073.43 

-69.26 

-1.02 

0.08 

12154.06 

-69.47 

-1.02 

0.09 

12381.12 

-70.07 

-1.04 

0.12 

10 

13270.68 

-77.22 

-1.33 

0.36 

13350.08 

-78.03 ■ 

-1.93 

-0.22 

13573.70 

-80.31 ■ 

-3.61 

-1.84 

11 

14453.74 

-85.31 

-1.68 

0.77 

14533.22 

-85.41 - 

-1.56 

0.92 

14757.10 

-85.66 ■ 

-1.21 

1.35 


excitations, we found it important to inclnde the experimentally derived vibrational energies reported by I Ghosh et al.l (119811 ) . 
although with fitting weights smaller by a factor about 10 than the rest of the data. Since DPotFit is not able to fit to energies 
directly, we used these energies to generate pseudo-transition frequencies, which allowed us to use data from higher vibrational 
states, up to n = 11, and obtain a much more reasonable fit at the higher energy part of the PEC. 

The values for the equilibr i um d istance re and dissociation energy Dg were fixed to the experimentally derived values of 
Gineeric 1 1969h: Cazzoli et al. ( 200fil) . After some experimentation the exponent o in Eo. m was set to 4, which allowed us 


to obtain an accurate fit to the higher energy levels. Four expansion variables Pi were used in the final results. Initial values 
for the parameters Pi were set to those obtaine d for the scaled PE G-S. It should be noted, however, that the final results are 
not very sensitive to the initial ab initio curve ( Barton et al. 2013h . The final empirical parameters are given in Table |4] 

We did not consider explicitly non-adiabatic corrections as these effects are expected to have a small effect (< 0.05 cm“^) 
on most transition frequences because of the relatively large nuclear masses of PN and the absence of low-lying excited 
electronic states. 

Table [3] gives residuals for a selection of energy levels. These results show that PEC-R generally improves on PEC-S in 
terms of agreement with experiment and all levels are reproduced to better than 1 cm~^. Table [5] gives residuals for very 
accurate pure rotational transitions in the microwave region for ®^P^^N and The ®^P^®N microwave data in this table 

are the only available experimental data for this isotopologue. Because of cancellation of errors calculated rotational transition 
frequences are for rotational transitions are much more accurate than for rovibrational ones. The ab initio PEC-A gives the 
worse results, with a root-mean-square (rms) error of 0.03 cm“^; PEC-S improves very considerably over PEC-A and gives 
an rms error of 0.001 cm“^. Finally, the empirical PEC-R gives the smallest rms of 7x10“® cm“^. 

Overall the empirical PEC-R reproduces the experimental transition wavenumbers used in the fit (excluding the w > 4 
simulated transitions) with an rms error of 0.004 cm“^. 

Fig. [T] compares the refined PEC-R with the ab initio PEC-A and the one-parameter semiempirical PEC-S. Due to the 
lack of experimental data for energies higher than about 14 500 cm“^, only the range r « 1.28 A to r « 1.80 A of the 
empirical PEC-R is well characterised. As one can see from Fig. [T] while within this range PEC-R and PEC-S agree to about 
10 cm“^ outside it the two PECs differ by up to about 400 cm“^. In particular, the peak visible for both PEC-A and PEC-S 
at r « 2.30 A corresponds to the minimum of the co curve (coefficient of the CASSCF wave function in the MRCI expansion, 
see discussion in Section I3.1|l . A more detailed study of the electronic structure of PN at long bond lengths would be needed 
to ascertain whether this feature has a physical origin or is an artifact of our calculation method, but this is beyond the scope 
of this work. We can however note that the difference in J = 0 energy levels between PEC-S and PEC-R for n > 11 grow 
approximately as 0.22 x (w — 11)^ cm“^ up to n = 47 (PEC-R energy levels are the higher ones) and reaches the maximum 
value of 210 cm“^ for that value of v. 

Despite this large difference in energy levels transition frequences with An = 1 (the strongest) and n > 11 show maximum 
deviations between PEC-R and PEC-S of only up to 10 cm“^. 
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6 Yorke et al 


Table 4. Parameters defining the empirical potential surface PEC-R: expansion coefficients /3i, dissociation energy De in cm equi¬ 
librium bond length Ve in A, and parameter p, see eq. For r < Ve the coefficients /5(3) and /5(4) are set to zero. 


Parameter 

Value 

Z)e/cm“l 

51940 

re/ /AA 

1.4908696082 

m 

2.2185248572 

m 

0.17749675548 

m 

0.16962236112 

m 

0.14242140368 

m 

0.22164327612 

p 

4 


Table 5. Differences of experimental llWvse et alJIlQ^ : ICazzoli et ani2006l) . (‘obs.’) pure rotational transitions with calculated ones 
(‘calc.’) for three potential energy curves PEC-X, X =A, S or R. All transitions have Au = 0 and J" = J' -j- 1. Energies are in cm“^. 


V 

J' 

obs. 

(obs. ■ 

- calc.) xlO® 




A 

S 

R 




pl4N 



0 

0 

1.567427 

-295.5 

10.9 

-0.5 

0 

1 

3.134828 

-590.9 

21.8 

-1.0 

0 

2 

4.702176 

-886.5 

32.6 

-1.6 

0 

3 

6.269446 

-1182.0 

43.5 

-2.2 

0 

4 

7.836611 

-1477.6 

54.3 

-2.7 

0 

5 

9.403643 

-1773.5 

65.0 

-3.5 

0 

6 

10.970523 

-2069.0 

76.1 

-3.8 

0 

7 

12.537217 

-2364.8 

86.9 

-4.4 

0 

10 

17.235939 

-3252.9 

119.4 

-6.0 

0 

11 

18.801641 

-3549.0 

130.4 

-6.4 

0 

12 

20.367025 

-3845.6 

141.1 

-7.2 

0 

13 

21.932071 

-4142.2 

151.9 

-7.7 

0 

14 

23.496750 

-4439.0 

162.7 

-8.3 

0 

15 

25.061036 

-4736.0 

173.3 

-9.0 

0 

16 

26.624905 

-5033.1 

184.2 

-9.5 

1 

1 

3.112628 

-606.1 

21.6 

-4.7 

1 

2 

4.668876 

-909.3 

32.3 

-7.1 

1 

3 

6.225046 

-1212.3 

43.2 

-9.4 

1 

4 

7.781110 

-1515.6 

53.9 

-11.8 

1 

5 

9.337046 

-1818.6 

64.9 

-13.9 

1 

6 

10.892821 

-2122.0 

75.6 

-16.4 

2 

2 

4.635485 

-932.7 

31.7 

-7.9 

2 

3 

6.180523 

-1243.8 

42.2 

-10.6 

2 

4 

7.725463 

-1554.2 

53.4 

-12.5 

3 

3 

6.135871 

-1276.5 

40.3 

-6.1 

3 

4 

7.669644 

-1595.4 

50.7 

-7.4 

4 

3 

6.091073 

-1310.7 

37.2 

2.6 




pl5N 



0 

1 

2.991601 

-561.6 

22.9 

1.1 

0 

2 

4.487339 

-842.8 

34.0 

1.4 

0 

7 

11.964490 

-2248.2 

90.6 

3.6 

0 

8 

13.459445 

-2529.4 

102.0 

4.2 

0 

10 

16.448686 

-3092.3 

124.7 

5.2 

0 

11 

17.942926 

-3373.8 

136.1 

5.7 

0 

12 

19.436878 

-3655.7 

147.2 

6.0 

0 

13 

20.930522 

-3937.6 

158.6 

6.5 

0 

14 

22.423832 

-4219.6 

169.8 

7.0 
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Figure 1. The thick, blue curve is the 8-parameter empirical potential PEC-R defined by Eq. and Table [4] (vertical scale given by the 
axis on the left hand side). The dashed, green curve is the difference between the ab initio PEC-A and PEC-R and the thin, red curve 
the difference between the one-parameter semi-empirical PEC-S and PEC-R (vertical scale given by the axis on the right hand side). The 
inset on the right is a blow-up of the area r = 1.28—1.80 A contained in the black rectangle, corresponding to the range characterised by 
the available experimental data; while within this range PEC-S and PEC-R agree within about 10 cm~^, outside this range differences 
can be as large as 400 cm“^. 



Figure 2. The ab initio MRCI-|-Q/aug-cc-pCV6Z dipole moment curve (DMC) of PN. 


3.4 Dipole moment curve (DMC) 

There appears to be no full DMC for the ground electronic state of PN available in the literature. Furthermore, there are no 
experimental absolute transition intensities, so it is difficult to judge the accuracy of a DMC. We therefore used the highest 
level of ab initio theory considered in this work, MRCI+Q / aug-cc-pCV6Z, to compute a new DMC. Dipoles were computed 
as the derivati ve of the MRCI+Q ene rgy with respect to an external electric field along the internuclear axis for vanishing 
field strength 1 Lodi Sz Tennyson 2Q10l) : we used field strengths ±5 x 10 a.u. and computed dipoles on a grid r = 2.1 to 
6.1 ao in steps of 0.05 ao. The corresponding DMC is shown in Fig. 0 This DMC has an eq uilibrium value of 2.739 D , whic h 
is close to the value of the dipole moment, 2.7465T0.0006 D, determined experimentally bv iRavmonda fc Klempereil lll97lll . 


3.5 Computing the line list 

The empirical PEC and the ab initio DMC described above were used to generate the line lists for and ®^P^®N with 

LEVEL 8.0. The integration range was chosen as r = [1.0,10.0] A and the grid comprised 45 000 points. We estimated the 
numerical error in computed energy levels to be less than 0.0005 cm~^. Table [6] summarises the results, showing the maximal 


© 0000 RAS, MNRAS 000, 000-000 























8 Yorke et al 


Table 6 . Summary of our line lists. 



31pl4N 

sipiSN 

Max V 

66 

68 

Max J 

357 

366 

Max V cm~^ 

51928.7 

51937.0 

number of lines 

692019 

743114 

number of energies 

13949 

14623 


Table 7. Sample extracts from the energy and transition files for The energy files contain 13 949 entries for and 14 623 

for while the transition files contain 692 019 entries for 3ipi4^ 743 334 fQj. 3lpl5^ 


N 

Energy file 

E g J 

V 

I 

F 

Transition file 

1 

0.000000 

4 

0 

0 

2 

1 

3.0148E-06 

2 

1.567437 

12 

1 

0 

3 

2 

2.8941E-05 

3 

4.702269 

20 

2 

0 

4 

3 

1.0465E-04 

4 

9,404461 

28 

3 

0 

5 

4 

2.5722E-04 

5 

15.673929 

36 

4 

0 

6 

5 

5.1376E-04 

6 

23.510567 

44 

5 

0 

7 

6 

9.0135E-04 

7 

32.914245 

52 

6 

0 

8 

7 

1.4470E-03 

8 

43.884806 

60 

7 

0 

9 

8 

2.1778E-03 

9 

56.422067 

68 

8 

0 

10 

9 

3.1207E-03 

10 

70.525819 

76 

9 

0 

11 

10 

4.3027E-03 

11 

86.195825 

84 

10 

0 

12 

11 

5.7506E-03 

12 

103.431825 

92 

11 

0 

13 

12 

7.4913E-03 

N: State counting number. 


I: 

Upper state counting number. 

E: State energy in cm 

1 


F: 

Lower state counting number. 


g\ State degeneracy. Aif: Einstein A coefficient in s 

J: State rotational quantum number. 
v: State vibrational quantum number. 


values of u, J as well as the wavenumber range considered in this work, which is far more extensive than any previous study. 
However, as discussed in Section our empirical potential curve is not well characterised for energies higher than about 
14 500 cm“^ and therefore transition frequences involving energy levels with energies higher than « 14 500 cm“^ may have 
large errors of up to a few hundreds of cm“^. 

Following Tennyson et al. 1 201,1 ) our line lists consist of an ‘Energy’ and a ‘Transition’ file, extracts from which are given 
in Table [3 The full line list can be downloaded from the Exomol web site www. exomol. com or from the VizieR service from 
the web site of the Strasbourg Astronomical Data Center cds.u-strasbg.fr, A small computer program which uses these 
files to compute spectra is given in the supplementary material. 


3.6 Partition Function 


The partition functions, Q{T), for and were calculated by direct summation over the energy levels using 


( 6 ) 


QiT) = 5ns (2 Ji + 1) ^ 

i=0 

where Ji is the rotational angular momentum quantum number and Ei is the energy of the of state i. Partition functions can 
be used to compute a variety of thermodynamic data. 

Exo Mol follows the HIT RAN convention of explicitly including the full atomic nuclear spin in the molecular partition 
function ( Fischer et al. 2003l 'l . In our case gns = 4 for and 5ns = 6 for 


The ^^ P^^N partitio n function is illustrated in Table [S] where it is compared to the values obtained from the parameters 


pr ovided by Rwir] ( 1981 1. The agreement is good especially considering a rather limited leve l of accuracy of the model used 
bv llrwinl ( 198ll l. At lower temperatures our values do deviate somewhat from those derived bv lirwirj ( 198l|), which is expected 
since the latter is only designed for T > 1000 K. In total we have considered temperatures from 10 K-5000 K in this work, and 
a full table of values for the partition function values on a fine grid for the two isotopologues can be found in the supplementary 
material. We also computed partition functions using energy levels obtained from the semi-empirical PEC-S to estimate the 
error in our computed Q{T). The conclusion is that our partition functions should have a relative error of less than 0.01%, 
i.e. they are correct to four significant digits. 
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Table 8. Tabulations of the partition function for 3lpl4pj given temperatures in kelvin compared to the values by lirwi 3 lll98lh . The 
latter were scaled by the gns = 4 factor. The full table can be found in the supplementary material. 


T 

This work Jrwin (T9M) 

100 

356.2 

515.6 

200 

711.2 

848.3 

300 

1068.2 

1156.7 

400 

1434.1 

1487.1 

500 

1817.9 

1848.6 

600 

2227.0 

2244.5 

700 

2666.1 

2676.0 

800 

3138.3 

3143.9 

900 

3645.3 

3648.7 

1000 

4188.4 

4190.8 

1500 

7468.0 

7470.7 

2000 

11717.8 

11719.8 

2500 

16964.4 

16961.2 

3000 

23229.7 

23215.1 

3500 

30535.8 

30500.1 

4000 

38905.7 

38834.1 

4500 

48363.9 

48234.7 

5000 

58936.7 

58719.8 



Figure 3. Overview of the absorption spectrum of 3lpl4yj g temperature of 1200 K. 


3.7 Examples of ro-vibrational spectra 

We used the new line l ist of to compute a number of absorption spectra in the form of integrated cross-sections as 

described by Hill et al. ( 20131 ). Figure [37fl shows the absorption spectrum of at T = 1 200 K. Only is presented 

as the difference wit h spectra for ^^P^^N can o nly be distinguished at high resolution. Figure 13.71 compares a spectrum of 
3ipi4N recorded by Ahmad fc Hamilton I lOOSh with the spectrum from this work. Due to pollution of the experimental 
spectrum with water lines, the spectra are far from identical; however they show similar structure, with the P and R branches 
clearly visible in both spectra. This demonstrates the accuracy of our results. Another illustration is presen ted in Fig. El where 
the microwave spectrum of collected in the Cologne Database for Molecular Spectroscopy (CDMS) ( Muller et al.llio^ ) 

compared to our spectrum. The agrement is excellent. 


4 DISCUSSION AND CONCLUSION 

We produced three new potential energy curve (PEC) for PN: an ab initio PEC-A, a one-parameter semi-empirical PEC-S 
and a fully empirical 8-parameter PEC-R. The produced PECs are applicable to both the and ®^P^®N isotopologues. 

The empirical PEC-R reproduces all known ^^P^^N experimental data up to u = 4 with a typical error of less than 0.01 cm“^ 
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wavenumber / cm’’ 

1280 1300 1320 1340 1360 


V 1.0x10^ 

0) 



0.96 - 


■ * I I I I I r—I I I 1 I , . I r*--T I 

1270. 1295. 1320. 1345. 1370. 

Wavenumbers (cm"') 

Figure 4. Comparison between the experimental spectrum of bv lAhmad &: HamiltonI l|l995l') and this work, at room temperature. 

Our work is given in the upper panel and has intensities in absolute units. The experimental d ata is given as a transmission spectrum 
and shows many strong features due to water. (Reprinted from Ref. llAhmad fc HamiltonlllOOfji . Copyright 1995, with permission from 
Elsevier.) 


and is expected to be the most accurate of the three. Transition frequencies for the minor isotopologue are expected 


to be somewhat less accurate. Experimental pure rotational transitions for both isotopologues with u < 4 (ICazzoli et al. 


20061 : Ahmad &: Hamilton 1995 ) in the microwave region are reproduced within 10 cm Levels with ii = 5tou = ll are 


expected to have larger errors, in the range 0.1 to 1.0 cm“ . The accuracy of semiempirical PEC-S should be in most cases 
not much worse than the quoted values for PEC-R. Errors in the ah initio PEC-A are larger by a factor 10-100. Levels with 
u > 11 are expected to have very large errors for all PECs, in the range 10-200 cm“^; however, transition frequences with 
II > 11 and An = 1 should be accurate to «5 cm“^. Further experimental and theoretical studies are needed to improve the 
accuracy for high v’s. 

An ab initio dipole moment was produced and, together with the PEC-R, used to compute line lists (comprising transition 
line positions and Einstein coefficients) for and Details of the line lists are reported in Table IH The line list 

cover all rotational-vibrational levels up to dissociation but, as discussed above, transitions involving u > 11 may be seriously 
in error. 

Partition functions a ccurat e to ab out 0.01% were also calculated for a range of temperatures for both isotopologues. 
Work done previously by lirwini (119811 ) shows a great deal of agreement with our partition function, with the exception of 
temperatures below 1000 K where Irwin’s results are not valid. 


5 SUPPLEMENTARY MATERIAL 

The Supplementary material to this paper includes the line lists and partition functions of ®^P^^N and as well 

as the PECs and DMCs used in these calculations together with the LEVEL 8.0 and DPotFit input/output files, and a 
Fortran program which uses the state and transition files to compute spectra. The line list are available from the Exomol 
web site www.exomol.com or from the VizieR service from the web site of the Strasbourg Astronomical Data Center via 
http://cdsarc.u-strasbg.fr/cgi-bin/VizieR?-source=J/MNRAS/xxx/yy 
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Figure 5. Compa rison between the room-temperature microwave spectrum from the Cologne Database for Molecular Spec¬ 

troscopy (CDMS) llMiiller et al.ll2005[l and this work. 


ACKNOWLEDGEMENTS 

This work is supported by ERC Advanced Investigator Project 267219. 


REFERENCES 

Ahmad I. K., Hamilton P. A., 1995, J. Mol. Spectrosc., 169, 286 
Bartlett R. J., Musial M., 2007, Rev. Mod. Phys., 79, 291 

Barton E. J., Chiu C., Golpayegani S., Yurchenko S. N., Tennyson J., Frohman D. J., Bernath P. F., 2014, MNRAS, 442, 
1821 

Barton E. J., Yurchenko S. N., Tennyson J., 2013, MNRAS, 434, 14691475 

Bernath P., 2005, Spectra of Atoms and Molecules, 2nd edn. Oxford University Press 

Cazzoli G., Cludi L., Puzzarini C., 2006, J. Molec. Struct. (THEOGHEM), 780-81, 260 

Goquart B., Prudhomme J. C., 1980, J. Phys. B: At. Mol. Opt. Phys., 13, 2251 

Goquart B., Prudhomme J. C., 1981, J. Mol. Spectrosc., 87, 75 

Gurry J., Herzberg L., Herzberg G., 1933a, Z. Phys., 86, 348 

Gurry J., Herzberg L., Herzberg G., 1933b, J. Ghem. Phys., 1, 749 

de Brouckere G., Feller D., Koot J. J. A., 1993, J. Phys. B: At. Mol. Opt. Phys., 26, 1915 

Dunning T. H., 1989, J. Ghem. Phys., 90, 1007 

Feller D., Peterson K. A., 2013, J. Ghem. Phys., 139, 084110 

Fischer J., Gamache R. R., Goldman A., Rothman L. S., Perrin A., 2003, J. Quant. Spectrosc. Radiat. Transf., 82, 401 

Ghosh S. N., Verma R. D., VanderLinde J., 1981, Can. J. Phys., 59, 1640 

Gingeric K. A., 1969, J. Phys. Ghem., 73, 2734 

Grein F., Kapur A., 1983, J. Mol. Spectrosc., 99, 25 

Hill C., Yurchenko S. N., Tennyson J., 2013, Icarus, 226, 1673 

Hoeft J., Tiemann E., Torring T., 1972, Z. Naturforsch. A, A 27, 703 

Irikura K. K., 2007, J. Phys. Chem. Ref. Data, 36, 389 

Irwin A. W., 1981, ApJS, 45, 621 

Kemeny A. E., Francisco J. S., Dixon D. A., Feller D., 2003, J. Chem. Phys., 118, 8290 
Le Floch A. C., Melen F., Dubois L, Bredohl H., 1996, J. Mol. Spectrosc., 176, 75 

Le Roy R. J., 2006, DPotFit 1.1 A Computer Program for Fitting Diatomic Molecule Spectral Data to Potential Energy 
Functions. University of Waterloo Chemical Physics Research Report CP-662R, http://leroy.uwaterloo.ca/programs/ 
Le Roy R. J., 2007, LE'VEL 8.0 A Computer Program for Solving the Radial Schrodinger Equation for Bound and Quasibound 
Levels. University of Waterloo Chemical Physics Research Report CP-663, http://leroy.uwaterloo.ca/programs/ 

Lodi L., Tennyson J., 2010, J. Phys. B: At. Mol. Opt. Phys., 43, 133001 
Lodi L., Tennyson J., Polyansky O. L., 2011, J. Chem. Phys., 135, 034113 
Lynas-Gray A. E., Miller S., Tennyson J., 1995, J. Mol. Spectrosc., 169, 458 


© 0000 RAS, MNRAS 000, 000-000 




































12 Yorke et al 


Maki A. G., Lovas F. J., 1981, J. Mol. Spectrosc., 85, 368 

McLean A. D., Liu B., Chandler G. S., 1992, J. Phys. Chem., 97, 8459 

Milam S. N., Halfen D. T., Tenenbaum E. D., Apponi A. J., Woolf N. J., Ziurys L. M., 2008, ApJ, 684, 618 

Miiller H. S. P., Schloder F., Stutzki J., Winnewisser G., 2005, J. Molec. Struct. (THEOCHEM), 742, 215 

Patrascu A. T., Tennyson J., Yurchenko S. N., 2014, MNRAS 

Peterson K. A., Dunning T. H., 2002, J. Chem. Phys., 117, 10548 

Pfeiffer B., Venkataramaniah K., Czok U., Scheidenberger, 2012, ADNDT, 100, 403 

Raymonda J., Klemperer W., 1971, J. Chem. Phys., 55, 232 

Saraswathy P., Krishnamurty G., 1987, Pramana, 29, 53 

Szalay P., Muller T., Gidofalvi G., Lischka H., Shepard R., 2012, Chem. Rev., 112, 108 
Tennyson J., 2014, J. Mol. Spectrosc., 296, 1 

Tennyson J., Hill C., Yurchenko S. N., 2013, in AIP Conference Proceedings, Vol. 1545, 6*^ international conference on 
atomic and molecular data and their applications ICAMDATA-2012, AIP, New York, pp. 186-195 
Tennyson J., Yurchenko S. N., 2012, MNRAS, 425, 21 
Turner B. E., Bally J., 1987, ApJ, 321, L75 

Turner B. E., Tsuji T., Bally J., Guelin M., Cernicharo J., 1990, ApJ, 365, 569 

Verma R. D., Ghosh S. N., Iqbal Z., 1987, J. Phys. B: At. Mol. Opt. Phys., 20, 3961 

Viana R. B., Pereira P. S. S., Macedo L. G. M., Pimentel A. S., 2009, Chem. Phys., 363, 49 

Visscher C., Lodders K., Fegley, Jr. B., 2006, A&A, 648, 1181 

Surkus A. A., Rakauskas R. J., Bolotin A. B., 1984, Chem. Phys. Lett., 105, 291 

Werner H.-J., Knowles P. J., Knizia G., Manby F. R., Schiitz M., 2012, WIREs Comput. Mol. Sci., 2, 242 

Wong M. W., Radom L., 1990, J. Phys. Chem., 94, 638 

Woon D. E., Dunning Jr. T. H., 1995, J. Chem. Phys., 103, 4572 

Wyse F., Gordy W., Manson E., 1972, J. Chem. Phys., 57, 1106 

Yadin B., Vaness T., Conti P., Hill C., Yurchenko S. N., Tennyson J., 2012, MNRAS, 425, 34 
Yamaguchi T. et al., 2011, Publ. Astron. Soc. Jpn., 63, L37 

Yamaguchi T. et al., 2012, Publ. Astron. Soc. Jpn., 64, 105 

Ziurys L. M., 1987, ApJ, 321, L81 


© 0000 RAS, MNRAS 000, 000-000 


potential energy / cm 



1 










